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Abstract 

This paper is concerned with modeling the dynamics of social outbursts of activity, such 
as protests or rioting activity. In this sequel to our work in [5], written in collaboration with 
J-P. Nadal, we model the effect of restriction of information and explore the effects that it has 
on the existence of upheaval waves. The systems involve the coupling of an explicit variable 
representing the intensity of rioting activity and an underlying (implicit) field of social tension. 
We prove the existence of global solutions to the Cauchy problem in as well as the existence 
of traveling wave solutions under certain parameter regimes. We furthermore explore the effects 
of heterogeneities in the environment with the help of numerical simulations, which leads to 
pulsating waves in certain cases. We analyze the effects of periodic domains as well as the 
barrier problem with the help of numerical simulations and discuss many open problems. 


1 Introduction 

The need to understand the spatio-temporal dynamics of social outbursts, such as protests or rioting 
activity, has been highlighted with many current events. For example, the fatal shooting of Michael 
Brown in Ferguson, Missouri (USA) by a white police officer in August of 2014 led to a two-week 
period of protests, which subdued until the grand jury’s decision not to indict the officer in question 
(Qj. A similar wave of protests was later observed in New York City (USA) after the grand jury’s 
decision not to indite the police officer involved in Eric Garner’s chokehold case m As many 
individuals throughout the U.S. and the rest of the world are keeping a close watch on these and 
similar cases there is a risk of increased social tension spreading throughout the country and abroad 
mmm- We think of protests as outbursts of social activity which are induced by an initial shock 
and grow, due to a self-reinforcement mechanism and/or continuing external shocks, for some period 
of time. We sometimes refer to this initial shock as the triggering event. Clearly, it is important 
not only to understand when these outbursts of activity will occur, but also how they will spread 
geographically. This motivated us to introduce various models that include what we believe to 
be the bare essential mechanisms which are necessary to capture the stylized facts of outbursts of 
social activity, in particular related to protests or civil unrest, in [5] written in collaboration with 
J-P. Nadal. The models we introduced involve the coupling of an explicit variable representing 
the intensity of rioting activity and an underlying (implicit) field of social tension. The effects of 
exogenous and endogenous factors as well as various propagation mechanisms are also included in 
the system. It must be said that the purpose of these series of works is not to understand the 
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economic, social or political origins of riots, and much less to discuss the legitimacy or lack thereof 
of any riot. These issues are much too complex to be studied with mathematical models which 
are obtained by making simplifying assumptions. Instead, our objectives are of an exploratory 
nature, where we wish to capture some qualitative facts of riots using our model with the aim of 
gaining insight into how riots spread and what the impact of restriction of information and various 
communications networks is on the spread and duration of rioting activity. 

The models of [5] introduce new and interesting mathematical structures for which fundamental 
questions, such as whether the systems have global solutions, have yet to be studied. The existence 
of traveling wave solutions in the case of non-local diffusion is also an open question, which we 
hope to address in future work. The purpose of this manuscript is four-fold. First, we introduce a 
model that captures the effect that restriction of information has on the success or failure of a riot. 
Second, we explore the existence of traveling wave solutions (upheaval waves), which provides insight 
into how fast riots spread. Third, we explore the effects of various heterogeneous environments on 
the success or failure of riots through numerical experiments. Finally, we take the first step in 
developing the mathematical theory called for by the family of models introduced in |5] , by proving 
the existence and uniqueness of solutions to the Cauchy problem. Furthermore, we discuss many 
open problems that are of interest from both the point of view of mathematics and the application. 
The approach of this work bears similarities to a recent literature on the use of mathematics in the 
analysis of uncivil and criminal activities [126, 123] . showing for instance that patterns, which are 

useful to understand, emerge when looking at the macroscopic scale. 

As the family of models we introduced in our work with Nadal [5] capture the effects of geographic 
proximity as well as social connections between different regions that are prone to protests or rioting 
activity, they are natural candidates to shed light onto answer to questions such as why some 
protest or even revolutions skip certain countries. There is evidence that restriction of information 
hinders revolutions and that increased social media access and technology has the opposite effect 
[12) [El OQ ES]. Naturally, increased access to technology has decreased the ability of any state 
to isolate their populations from the rest of the world. The models we introduce here modify 
those from [5] in order to take into account the effect of restriction of information in a system. 
Under some parameter regime our model is able to capture “upheaval waves,” a disturbance of 
rioting activity or protests that propagates through space. These disturbances are sparked by 
a triggering event of sufficient intensity. It is of interest to quantify the speed at which these 
upheaval waves move, something that is facilitated by the models we introduce. Moreover, as the 
restriction of information is reduced we observe a shift in the qualitative behavior of the upheaval 
waves. Particularly, when there is some restriction of information there exists upheaval waves, 
under certain parameter regimes, which travel at a unique speed. However, as the restriction of 
information is lifted completely the waves can travel with speed c for c > c*. Presumably, the 
intensity and range of the initial shock will have an effect on the spreed. In fact, from numerical 
simulations we observe that initial levels of rioting activity of the form u(x , 0) = 0(e ~ kx ) lead to 
traveling wave solutions whose speed depends on k. This is analogous to what is observed in the 
well-known Fisher-KPP equation, see for example [23| and references within. 

There have been other models of rioting activity introduced in the literature, see for example [8] 
EllEj. We refer the interested reader to [5] for a summary of these works. One of our contribution 
is the introduction of a continuous model, which is amenable toward more rigorous mathematical 
analysis. In fact, our model takes a new direction aiming at using systems of PDEs, an approach 
which has been useful in gaining insight into other social phenomena, such as urban crime mm, 
social segregation [21], opinion dynamics [23j . 

Outline of the paper: In section [2] we introduce the basic model along with some variations which 
we study in this work. In section [3] we prove some preliminary results needed for future sections. 
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The existence of global solutions to the Cauchy problem in for d > 1 is discussed in section [4| 
In section [5] we discuss the existence of traveling waves solutions and conclude with a discussion in 
section 0 

2 Dynamics of social outbursts of rioting activity 

2.1 The basic model 

In this section we summarize the models of |5j. The derivation began with a network of nodes which 
represent centers which are prone to rioting activity. We approximate the global behavior of such 
network by deriving the appropriate continuum limit, which then enables us to perform rigorous 
mathematical analysis. Let HcK 2 represent the domain of interest, i.e. a city or country. For a 
location x E II and time t we let u(x, t ) represent the level of rioting activity and v(x, t ) represent the 
social tension in the system. The fundamental assumption in our model is a dichotomy of external 
influences and internal influences in rioting activity: an external shock is necessary to initiate 
the process (exogenous factor), but the outburst is maintained by a self-reinforcement mechanism 
(endogenous factors). As these outburst are short-lived one must also include mechanisms that lead 
to an eventual self-relaxation of the activity. For simplicity we summarize the assumptions that we 
make below and refer the reader to [5] for a more detailed derivation. 

1. Exogenous factors : there are external factors that can either trigger or rekindle an outburst 
of protests or rioting activity. An example of this would be the shooting of Michael Brown 
in Ferguson, Missouri (USA) or the acquittance of the police officers involved in the beating 
of Rodney King in 2001 in Los Angeles, California (USA). Such events are modeled as point 
sources that are placed at the time and location of the triggering events. For example, if n 
external events or shocks occur at places and times (sj,U) 6 HxR then the source term for 
the social tension is given by: 


n 

s(x, t) .— Aj ^ ( §t=ti,s=Si j 

i= 1 

where A; is the intensity of the i th event. For the sake of simplicity, we will always assume 
that Ai = Ao for all i = 1 ...n in what follows. 

2. Endogenous factors: once an outburst of activity is initiated, it is maintained by a self¬ 
reinforcement mechanism, up to a certain carrying capacity. This is analogous to the repeat 
and near-repeat victimization effect assumed in the models of Short et al. in |25j and Beresty- 
cki and Nadal in [4J- In essence, this is the phenomena that criminal activities in certain 
location lead to an increased probability of that same location falling victim to more crimi¬ 
nal activities. The self-reinforcement mechanism is modeled by the function G(z), which is 
included as a source term in the dynamics of the level of rioting activity and satisfies: 

G(z) > 0 for (0, zq). G( 0) = 0 and G{z) < 0 for z > zq. 

An example of this is a KPP-type term: G(z) = z( 1 — ^) for z € (0,zo) for some zo > 0, 
where zq is the carrying capacity. The need for a carrying capacity arises from the fact that 
there are a finite number of people that will protest or in the case of riots a finite number of 
cars and buildings that can be destroyed. 
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3. Critical threshold for the social tension: an outburst of rioting activity or protests only occurs 
if the social tension in the system is sufficiently high. We introduce a critical threshold pa¬ 
rameter a G R + which determines when the endogenous factors begin to play a role. In other 
words, when the social tension is above this critical value the self-reinforcement mechanism 
is initiated. This gives rise to a transition function that depends on the social tension r(u). 
This transition function can be, for example, a sigmoid function: 

, v 7 

r(z) = - W, -r, 

V ' 1 + e -dC-a) 

where (3 > 0 provides a measure of the transition between a system without self-reinforcement 
and a self-reinforcing system and 7 gives us the units of the inverse of time. In other words, 
it provides a measure of how fast the transition is between a system that does not include 
the endogenous factors and a system with the full-force of these factors. Note that large f3 
provides fast transitions. The coefficient a is the critical tension threshold, which provides the 
level of tension necessary for the self-reinforcement mechanism to begin. 

4. Self-relaxation: there is a natural decay rate for the level of rioting activity, denoted by the 
parameter c 0 . Moreover, the level of rioting activity presumably will have an effect on the 
social tension and so it is natural to assume that the social tension deceases at a rate which is 
proportional to the level of rioting activity. Specifically, the decay is given by —h(u)v where 

h(u) = — -r— 

(1 + mu)P, 

for some 6,p,m > 0. The parameter 9 measures the decay in the absence of any riots, p 
controls the effect that the level of rioting activity has on the social tension, and m provides 
the units of the inverse of the level of rioting activity. Note that when u is large h(u) is small, 
modeling the fact that high levels of rioting activity will lead to a slow decay of the social 
tension level. 

5. Geographic proximity: while geographic proximity does not provide the full picture, especially 
nowadays, it still play an important role. For example, if we assume that the level of rioting 
activity u(x,t ) is simply given by the fraction of people protesting then local diffusion is 
appropriate in modeling the spread of the activity. Geographic proximity can also have an 
effect on the spread of social tension. Indeed, one might feel the effects of what occurs in 
neighboring locations more significantly. This effect is modeled by the Laplace operator and 
we will denote by D \, D 2 > 0 the diffusivity coefficients of u and v respectively. 

6 . Social proximity: one can argue the social connection between two places might be even 

more important than geographic proximity. Two locations that are socially connected might 
be more likely to influence each other than two places that are neighbors. The importance 
of this consideration is seen, for example, in the fact that urban centers are more likely to 
influence smaller cities or towns than vise-versa. Thus, it is important to include these social 
connections. These are modeled though the interaction kernel, which provides a 

measure of how much influence location y has on location x. In a sense, fl(x,-) gives the 
domain of influence of x, that is, the strength of the influence that other locations have on 
x. On the other hand, J(-,y) gives the range of influence of location y, that is, the influence 
that location y has on all other locations. 
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Combining the above assumptions we state the most general model: 

ut(x, t ) = D\Au(x, t) + r(v(x, t))G(u(x , t)) — wu( x, t), (la) 

vt(x, t) = D 2 Av(x, t) + k J(x, y)v(y) dy - ( h(u(x , t)) - rj)v(x, t ) + s(x, t) + v b , (lb) 

Ja 

for x E 11 and t > 0, where fi C or ff = M. d . In (jTb]) v b is the base social tension value. Note that 
we model the level of rioting activity that exceeds the base rate. This will be accompanied with the 
initial conditions: 


u(x,0) = uq(x) and v(x, 0 ) = vg(x), 


( 2 ) 


where u q(x),vq(x) G L°°(fl) n L 1 (H) and no-flux boundary conditions if 11 is bounded and limiting 
conditions if 11 = That is, from now on if 11 is bounded we assume that 


t —u(x,t) = -—v(x,t) = 0 for x G dQ and t > 0, 
on on 


(3) 


where n is is the outward unit normal. Some of the theory, such as the existence of solutions to 
the Cauchy problem, will cover the cases when d > 1. The theory of traveling wave solutions will 
be restricted to d = 1 since our main interest is the study of planar traveling waves, which can be 
reduced to the understanding of the evolution of the wave profile in one dimension. 

Remark 1. Generally speaking the dispersal of social tension in system one (JT]) will be either strictly 
local or strictly non-local, that is D 2 > 0 and n = 0 or D 2 = 0 and k > 0. However, we state the 
system in its generality to provide a general global well-posendess theory for the Cauchy problem. 


2.2 Modeling restriction of information 

In order to understand the effect that restriction of information has on the ignition and spread of 
protests we modify system 0 to include the fact that a revolution will begin only if the number 
of protests are above a certain threshold. The paper m used this concept with the aim of finding 
a model for the dynamics of revolutions, such as the Arab Spring. The idea is that small protests 
will simply go unnoticed. Moreover, it is also reasonable to assume that the tighter the restriction 
of information the larger the protests have to be in order for them to gain any momentum. We 
choose to model the restriction of information with the reaction term G a (u) which is assumed to 
be of ignition-type: 

G a {z) = 0 for z G [0, or], G a ( 0) > 0 for z G (a, 20 ) and G Q (z) < 0 otherwise, (4) 

for a G [0, zo]. We introduce a to represent the level of restriction of information in a system with 
larger values of a corresponding to higher levels of restriction of information. Naturally, there are 
many questions of interest which can be explored, such as what the effect of technology reaching the 
hands of more and more people is on the existence and speed of upheaval waves. Mathematically, 
this corresponds to determining the effect of taking the limit a —> 0. Another question that can be 
explored is that of why revolutionary waves skip certain countries. Needless to say, our work can 
only shed some light on this phenomenon and we do not claim to provide definite answers. 


2.3 Non-local diffusion 

When the social connection between two locations is only a function of the distance between those 
locations then the interaction kernel has the form J{x,y) = J(x — y) and the system will have the 
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following form: 

ut(x, t) = D\Au(x, t ) + r(v(x, t))G a (v(x, t)) — uju(x, t ), (5a) 

v t (x, t ) = D 2 Av(x, t) + k / J(x — y)v(y) dy - (. h(u(x , f)) - f) + s(.t, f) + v b . (5b) 

Jn 

On the one hand, this restricts the types of social interactions that we consider significantly, on 
the other hand, system ([5]) is a useful simpler model that allows for the exploration of non-local 
influences. For example, one can determine the role that the range of influence plays on the initi¬ 
ation, duration, and spread of rioting activity. For this special class of kernels, it is reasonable to 
expect that system ([5j will exhibit traveling wave solutions, whereas there is no reason to expect 
such solutions for a more general class of kernels. 


3 Preliminary analysis 

3.1 Non-dimensional system 

Before analyzing ([I]) let us first write the dimensionless model in order to remove superfluous 
parameters. Now, as discussed in [5] it is natural to make the assumption that the timescale over 
which the endogenous effects are observed is smaller than the exogenous effect, which then implies 
that uj > 8. Under this assumption the natural time scale is 1 /uj. The characteristic length scale 
is y/D[ /w, which is roughly the distance over which the level of rioting activity diffuses during the 
characteristic timescale. A natural change of variables then is: 

I UJ „u UJ 

t = cut, x = . / —x, u = —, and v = — v. 

V D \ zo v b 

This gives the non-dimensional system: 

ut = Au + r(v)Ga(u) — u, ( 6 a) 

v t = DAv + k J(x, y)v(y) dy - ( h{u) - k 2 )v + 1 + s(x, t ), ( 6 b) 


where the transition, decay, and source functions are now given by: 


r ( z ) = , , _-*(*-* y h ( z ) = 


(1 + mz)P ’ 


and s(x, t) = A ^ 5 t =t ir 


=ti,X=Xi 5 


i—1 


1 + e“AO 

the non-dimensional parameters are: 

7 , Kv b r] ~ A 0 _ au /3v b D 2 _ _ 

p = —, k = -, k 2 = —,A = —, a = —, pi =- = —, m = mz 0 , a = a zq. 

UJ UJ UJ V b Vq UJ 


and Ga (it) = 0 for u 6 [0, a] and Ga{u) = (u — ci)(l — u) for u ^ [0, a] with a € [0,1]. For 
notational simplicity we replace a with a for the remainder of the paper. Note also that since 
G a was normalized we gained a dimensionless parameter p that measures the intensity of the self¬ 
reinforcement mechanism. It will soon become apparent that this parameter, which we refer to as 
the self-reinforcement parameter because it measures the strength of the endogenous factors, plays 
a significant role in the behavior of our model. 
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3.2 Constant steady-state classification 

In this section we classify the constant steady-state solutions of © in the case of local diffusion, 
i.e. k = 0. We rewrite the system for convenience as follows: 


ut = A u + <J>(u, v), 


(7a) 


vt = DAv + k 


J(x,y)v(y) dy + ^(u,v) + s(x,t), 


(7b) 


where, 

<h(u,u) := r{v)G a (u) — u and ^{u, v) := — (h(u) — k 2 )v + 1. 
Any constant steady-state solution of ([7]) with k = 0 must satisfy: 

1 


v*(u) = 


h(u) - k 2 

Note that v*(u) is an increasing function of u such that 

lim. v*(u) = oo, 


and <L(u*(-u), u) = 0. 


( 8 ) 


for u = -=- 


\U r ,Vr 


it 


7)‘ /r -i 


Note that u > 0 as k 2 < 1 and that for any constant steady-state solution 


rolds that u c < min{l,tz} and u*(0) < v c < u*(l). In order to be able to analyze 
the effect of the transition parameter /3 we normalize the critical social tension value a to be the 
average between u*(0) and u*(l). Then the number of constant steady-states will depend on the 
restriction of information parameter a and the interplay between the transition parameter /3 and 
the self-reinforcement parameter p. We describe the different regimes of the parameters below and 
illustrate them in Figures [I] and [2j The various characterizations of the regions in the figure were 
determined through a combination of analytical results and heuristics, and confirmed numerically. 
See the Appendix (section [8]) for more details. It is important to keep in mind that this analysis is 
for the case when r(u), h{u ) and G a (u ) have the specific forms discussed above. 


1. No restriction of information: for a = 0 there can be either one, two, or three constant steady- 
state solutions. As can be observed in Figure [l] there exist two curves, (p) and /3 2 (p), that 
separate the number of steady-state solutions. A third curve f3s(p) separates the system from 
being a monostable system, where there is only one stable constant steady-state solution, and 
a bistable system, where there are two stable constant steady-state solutions. Finally, the 
curve /? 4 (p) divides the region where there are three constant steady-state solutions into two 
regions. In one region the non-excited steady-state solution is globally stable and in the other 
region the fully-excited steady-state solution is globally stable. We provide a more detailed 
summary of all of the regions below. 

• Figure [l] Region I: This region corresponds to a weak self-reinforcement effect. That is, 
if p is small then the only constant steady-state is the non-excited state (0, u*(0)). All 
solutions of the system ([7]) approach (0,u*(0)) in the long-term, see Lemma [2| 

• Figure [l] Region II: This region corresponds to intermediate and strong reinforcement 
effect but slow to intermediate transitions. As p increases the existence of two or three 
constant steady-state solutions depends on the transition parameter (3. Slow transitions 
lead to the existence of two constant steady-state solutions. For /3 below fisip) the 
system can be classified as monostable with one stable steady-state solution and the 
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other unstable (this is region lib of Figure [l]). Transitions in an intermediate region 
(region Ha of Figure [T]) lead to a bistable system with two stable constant steady-state 
solutions. 

• Figure[l] Region III: This region corresponds to an intermediate to strong self-reinforcement 
effect and fast transition rates. For p sufficiently large system 0 has three constant 
steady-state solutions. These correspond to the non-excited state (0,w*(0)), the warm 
state , by which we mean a steady-state with a non-trivial amount of social tension and 
rioting activity but not the maximum level, and the excited state (u*,v*(u*)). The warm 
state is unstable in this regime and the non-excited and excited states are stable. For 
intermediate values of p (region Ilia) the non-excited state is globally stable. On the 
other hand, for p to the right of f5i(p) the fully-excited state is globally stable. 

2. Moderate restriction of information: for 0 < a < 1 there are generally either one or three 
constant steady-state solutions. The existence of two steady-state solutions does occur when 
the parameters lie on the curve /3i(p), which is illustrated of Figure [2j However, in this cases 
the second non-zero steady-state solution is degenerate as <F(u, v*(u), ) < 0 for all u > 0. The 
curve /?2 (p) separates the regions where F(u ) < 0 for all 0 < u < min {1, it} and F(u) > 0 for 
some u G (0, min{1, It}). 

• Figure [2] Region I: This region corresponds to weak self-reinforcement effects. This if, if 
p is small then the only constant steady-state solution is the non-excited state ( 0 , u*( 0 )). 
As in the previous case of weak reinforcement mechanisms we expect that all solutions 
of system Q approach (0, u*(0)) in the long-term, see Lemma [2j 

• Figure [2] Region Ila: This region corresponds to medium-range reinforcement effects. 
As p increases there is a bifurcation into a regime where three constant steady-state 
solutions exist. In this case the non-excited state and the fully-excited state are both 
stable; however, the non-excited state will be globally stable. 

• Figure [ 2 ] Region lib: This region corresponds to large reinforcement effects. As p increases 
further there is a bifurcation into a regime where the fully-excited state is now globally 
stable. In this case the non-excited state is the preferred state. 

3. Maximum level of restriction of information: when a = 1 then (0, u*(0)) is the only constant 
steady-state solution. 



Figure 1: Bifurcation diagram for p and fi of © when k = 0 and there is no restriction of information. 





Figure 2: Bifurcation diagram for p and /3 of ([d]) when k = 0 and there 0 < a < 1. 


3.3 Comparison principle 

Many of the results in this work will rely on the construction of suitable super-solutions and sub¬ 
solutions. In order to be able to use these solutions as barriers we need the following simple 
comparison principle. 

Lemma 1 (Comparison principle). Let and (rt 2 ,?; 2 ) be two solutions of 0 with initial data 

satisfying ni(x,0) > ^(agO) and ui(.x,0) > u 2 (;r,0) in for d > 1. Then ui(x,t) > w 2 (.T,f) and 
v\(x, t ) > u 2 (x, t ) for all t > 0. 

Proof. Let u = u\ — U 2 and v = v\ — V 2 , then (u, v ) satisfy: 

u t = Au + (r(vi)G' a (u) -l)u + r'(v)G a fu 2 )v, (9a) 

v t = DAv + k j J(x, y)v{y) dy + (fc 2 - h(u\))v - h'(u)v 2 U, (9b) 

for u, u, v > 0. Since r'{y)G a {u 2 ) > 0 and h'{u)v 2 < 0 the system is monotone and thus by standard 
theory we conclude that u(x,t),v(x,t) > 0 for all t > 0 . J3 

3.4 Liouville type results 

As the ultimate distributions of the level of rioting activity and the social tension in the system 
are important to understand, the question of whether the only steady-state solutions to 0 are the 
constant solutions is not only of mathematical interest but also of general interest. It is certainly the 
case that for parameters in regions I of Figure [I] and Figure [2] where the reaction terms are strictly 
negative, the only steady-state solutions will be (0,u*(0)). Moreover, in the case when the social 
tension does not diffuse we can prove that the only steady-state solutions are indeed the constant 
ones. 

Proposition 1 . Let D = k = a = 0 and p , /3 be chosen from region Ha or lib of Figure [IJ Let ( u , v ) 
be a classical non-negative bounded steady-state solutions to 0 with no-flux boundary conditions 
0 on a bounded domain then either (u,v) = ( 0 , u*( 0 )) or (u,v) = (u *, v*(u *)). 
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Proof. Under the conditions stated in the proposition any steady-state solution will satisfy: 



Define f(u) = r &a(u) ~ u for which, based on the chosen parameters, there exists a 

unique u* such that: 

/(0) = f(u*) = 0, f(u) > 0 for u G (0, u*) and f(u) < 0 otherwise. 

Without loss of generality assume let us assume that u ^ 0. By the strong maximum principle we 
conclude that u(x) > 0 for all igO. Let m = min^gQ u(x) > 0 and M = max Ig jj u(x) > 0. Since 
f(u) is positive between (0,rt*) then by the strong maximum principle and Hopf lemma we know 
that m > u*. On the other hand, f(u) < 0 for u > u* so we can conclude, in a similar fashion, that 
M < u* , which implies that u = u* and thus v = v(u*). ;Q 

It would be interesting to find the parameter regimes, if such exists, where the only steady-state 
solutions are the constant ones for the fully-parabolic system. For example, when the parameters 
lie in the region Ha and lib of Figure [lj so that system ([T]) with k = 0 is monostable, are the only 
steady-state solutions the zero steady-state, (0, u*(0)), and the excited steady-state, (u*,v*(u*)). 
This would be a natural extension from results that are known to hold for the single parabolic 
Fisher-KPP equation. 


4 The Cauchy problem 

This section is devoted to the study of the Cauchy problem in M. d for d > 1 with one triggering 
event. We will take the triggering event to be part of the initial condition. More specifically, we 
consider the system: 

ut = A u + r(v)G a {u) — it, for x G M d , t > 0 , ( 10 a) 

vt = DAv + k j J(x, y)v(y) dy — ( h(u ) — k 2 )v + 1, for x € M d , t > 0, (10b) 

u(x,0)=uq(x) and v(x, 0 ) = vq(x) + A5 x =x, for x G M d , t = 0 . ( 10 c) 

Let C;,(M rf ) denote the space of continuous bounded functions on M. d and consider the Banach 
space U = Cb(M rf ) x C),(M d ) endowed with with the L°°— norm. We obtain the following global 
existence result. 

Theorem 1 (Global solutions to the Cauchy problem). Let (uq(x),v o(x)) G U, D > 0, and 

sup / J(x,y)dy< oo. (11) 

ISM 1 * J«. d 


Then, there exists a unique, positive, and global solution in time (u(x,t),v(x,t)) €ld to (10). 


Proof. The proof is done in two steps. In the first step we prove the local existence of a solution in 
a suitable subspace of U. The second step is a simple continuation argument. 
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Step 1: (Local existence) The solution to (10) with D > 0 can be re-written as: 


u(x,t) = e t JC(x, t) * uq(x) + / e t+s K.{x, t — s) * r(v(x, s))G a {u(x, s)) ds, (12a) 

Jo 

v(x, t ) = JCd(x, t ) * uo(x) + / JCd(x, t — s)* Ik / J(x, y)v(y, s ) dy — (. h(u(x, s)) - k 2 )v(x , s) + 1 


/o 


(12b) 


constant D. If D = 0 then (12b) instead has the form: 


where 1C is the heat kernel with diffusivity constant one and /Co is the heat kernel with diffusivity 

(13) 


u(x, t) = v 0 (x) + J (^k J J{x,y)v(y,s) dy — (h(u(x,s)) — k 2 )v(x,s) + 1^ ds. 
Let T[(u,v)\ := (Ti(u, v), T 2 (u, v)) be a map on U defined by: 

T\(u, v ) := e~ t K.(x, t) * uq(x) + / e~ t+s IC(x, t — s) * r(v(x, s))G a (u(x, s)) ds, 


(14) 


/ o 


T 2 {u,v) := ICd{x, t) *vo(x) + / K,D{x,t — s)ds 

Jo 

ft 


+ J Ko{x,t - s) * J J{x,y)v(y,s) dy - (h(u(x,s)) - k 2 )v(x,s)^j ds. (15) 


If D = 0 then we instead set: 

T 2 {u,v) := v 0 (x) + J [k j J(x,y)v(y, s) dy - (h(u(x, s) - k 2 )v(x, s) + 1 ) ds. 


(16) 


Our goal is to find a subset of U that is invariant under the transformation T. For this purpose, 
given T, R> 0 define the set: 


Xr,t = { (u, v) E U : || (u, v) || Xt < i?| , 


where, 


\\{u,v)\\ Xt = ||u(i)-e t K.(x,t)*u 0 {x)\\ oo + \\v(t)-v 0 {x)-t\\ oo , 
if D = 0. On the other hand, if D > 0 then we instead use the norm: 


II(^)ILy t = \\u{t)-e t JC(x,t)*u 0 (x)\\ oo + 


v(t) — ICd(x, t) * vq(x) — / Kd(x, t — s) ds 


Take (u, v ) E Xr t T, to show invariance we obtain the following estimate with the use of Young’s 
inequality for convolutions: 


\Ti(u,v)-e t X(x,t) * ti 0 (x)|| < [ e t+s ||/C(x, t - s) * r(v(x, s))G a (u(x, 

Jo 


s ))lloc ds 


< e t+s \\JC(x,t - s)!^ ||r(u(x,s))G Q ,(rt(x,s) 
Jo 

< C [ e~ t+s \\r(v(x, s))G a (u(x, s))^ ds 

Jo 

< CR(1 - e-% 


ds 


( 17 ) 


ds. 
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as || r(v(x, s))G a (u(x, s))|| < \G' a \ ||tt|| < CR. Similarly, for for D > 0 we obtain the estimate: 


(it, v) - JCd{x, t) * vq(x) - / ICd{x, t — s) ds 

Jo 


< C 


J{x,y)v(y,t)dy 


ds 


— / ( h(u(x,s )) — k 2 )v(x,s) ds 

Jo 

< C sup / J(x,y) dy t 

\xeIR d J M d / 

< Gift, (18) 

where we have used condition © and the fact that for ||'ifc|| 00 sufficiently small then h(u) — k 2 > 0. 
From © and ( |18[ ) we observe that for T sufficiently small T[u, v] G Xt.r- To show that T is a 
strictly contractive map, take (ui,ui), ( U 2 ,V 2 ) G Xt,r, we obtain 

||Ti(ui,ui) - 71 (« 2 ,^IL < |G' q | sup ||ui-u 2 || 00 (l-e“*), 

0 <t<T 


||75(iti,ui) - T 2 (u 2 ,v 2 ) 11^ < sup / J(x,y) dy + C] sup ||ui - t + Csup Wu! - u 2 \\ 00 t. 

ya:eR d 7ffi d J 0<t<T 0<t<T 

Therefore, for T sufficiently small there exists and 0 < k < 1 such that: 

||T[ui,ui] - T[u 2 ,v 2 }\\ u < k ||(wi,ui) - (u 2 ,v 2 )\\ u ■ 

Given that T is a strict contraction then there exists unique (u(x,t),v(x,t)) G U, which is a fixed 
point of T and thus satisfies ©• 

Step 2: (Continuation of the solution) Recall that we are looking for positive solutions. A simple 
continuation argument allows us to conclude that either T = oo or 

lim u(x,t) = oo or lim v(x,t) = oo. (19) 

t-*T~ t->T~ 

Since, ut < 0 for any u > u* we see that ||rt(x, i)^ < max{||rio(x )|| 00 ,u*}. On the other hand, 
v(x,t) < vo(x)e Mt where 


M := sup / J(x,y) dy. 


This rules out (19) and thus the solutions exists for all time. 


a. 


5 Traveling wave solutions without non-local effects 

From the point of view of the modeling it is of particular interest to look at the invasion of protests 
into regions that are experiencing low levels of protesting or rioting activity. A benefit of modeling 
riots with system ([b]) is that a wave of rioting activity or protests that spread geographically is 
represented by traveling wave solutions , i.e. solutions that propagate in a given direction with 
a constant profile. In the context of the model we refer to these solutions as upheaval waves or 
protests waves. In fact, for the local model (k = 0) we can determine what parameter regimes lead 


to such solutions for system (10) with k = 0. To be more precise, an upheaval wave is given by 


u(x, t ) = 4>(x — ct ) and v = ip(x — ct ) with c G M which satisfy: 

<j>"(z) + cxjJ(z) + ^(z),^(z)) = 0, 

+ cip'{z) + ’$>((j)(z),'ip(z)) = 0, 

0 < <j)(z) < u*, u* ( 0 ) < ip(z) < 

4>( — oo) = it*, 0(+oo) = 0 ,t/>( — oo) = V*(u*), tj}(+oo) = v*(0), 


( 20 ) 
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for z = x — ct. In system (20), (u*,v*(u*)) represents the excited constant steady-state solution and 
(0, v*(0)) represents the non-excited steady-state solution, with no rioting activity and the lowest 
social tension level. It is useful to define the following function: 


9 “ W ;=r (kW^) c “ W - 

as the speed of the upheaval waves will depend on the quantity: 

F(u) := f g a {s) ds. (21) 

Jo 

We first state a result about the eventual extinction of rioting activity when the parameters are 
chosen from regions I of Figure[l]and Figure[2j Of course, this implies the non-existence of upheaval 
waves. 


Lemma 2 (Extinction of rioting activity). Let k = 0, a = 1 or 0 < a and p , (3 chosen from regions 
I of Figures [7]{^. For any initial data (uq(x), vq(x)) £ L 1 (M. d ) x L 1 (M. d ) and source term: 

N 

s(x,t) = a y ] $t=ti,x=xii 
1=1 

the solution (u(x,t),v(x,t)) to Q approaches (0, u*(0)) as t oo. Specifically, 

u(x, t) —> 0 and v(x, t) —> u*(0), 


as t —>• oo, uniformly in x. 

Lemma [2] has the simple interpretation that in a fully censored state upheaval waves will never 
exist. At the same time even as the restriction of information decreases, if p is not sufficiently large 
there is also no possibility of an upheaval wave taking off. How large p has to be depends on the 
restriction of information; the higher the restriction of information the higher p needs to be in order 
for a riot to get momentum. 

Proof. (Lemma [2]) First assume that there is only one triggering event, in which case we assume 
that the source term can be taken as part of the initial condition and s(x,t ) = 0. Let (3 and p 
satisfy the hypothesis of the lemma, there exists an r > 0, which is sufficiently small, such that 
r(v)G a [u ) — u < —tu. For example, when a = 1 then r = 1. Let 

u(x, t ) = e~ Tt uo(x ) * KL(x, t), ( 22 ) 


where K is the heat kernel. Note that (22) satisfies ut = A u — tu. Thus, any solution u(x,t ) of 0 
is bounded above by u(x, t ) for all t > 0. Moreover, we have the bound 


u i x i t) lloo — e T t|«o(x)*/C(x,i)|| 00 <C7(||«o(®)|| 00 ,||X:|| 1 )e 


—rt 


In particular, ||u|| — > 0 as t — > oo uniformly in x. Furthermore, note that h(u) > /i(||w|| ) since 

h(u) is a decreasing function of u. Let w = v — v*(0) = v — h *^_ ko , then w satisfies: 


wt = D Aw — (h{u ) 


k 2 )w + 


h(u) - k 2 \ 

Mo) - fo ) 


+ s(x,t). 
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Using (22) we obtain that 


lim 

t—* OO 


1 - 


h(u) - k 2 


= 0 , 


h{0 ) - k 2 

also uniformly in x since u converges uniformly. Thus, we observe that ||w;(:e, t)\\ —> 0 as t —> oo 
uniformly in x. The result for the case when there are N triggering events follows by the argument 
above by taking the initial conditions to be the uq{x ) = u(x,tjy) and vq(x) = v(x,tw). □ 

On the other end of a the spectrum, as the restriction of information begins to be lifted there are 
certain parameter regimes that will allow the existence of revolutionary waves. This is expressed in 
the following result. 

Theorem 2 (Upheaval waves and restriction of information). Let 0 < a. < 1. 

1. For p and f3 in regions Ha, Ilia, or Illb of Figure [7] or regions I la and lib of Figurethe 
system (20) admits a unique solution ( 4>(z),'if(z),c *) up to translations. Furthermore, 


(i) sgn{c*) = sgn(F(u*)), where F{u ) is defined by (21). 

(ii) <j/(z) < 0 and ip'(z) < 0 for all \z\ < oo. 


2. If a = 0 and p and (3 lie in region lib of Figure then ( 20 ) admits a unique c* G M + and 
solutions (cf)c(z) , if c {z)) with z = x — ct for any c> c*. On the other hand, when c < c*, such 
solutions do not exist. 

In contrast to Lemma [2j Theorem [2] tells us that as the restriction of information decreases in 
a state the possibility of the existence of upheaval waves begins to increase. In fact, when a = 0 
there exists a wide parameter regime in the p — /3 plane that allows for the existence of traveling 
waves. Even more interesting is the possibility of upheaval waves that move with an arbitrarily 
large speed (depending on the initial shock). For 0 < a < 1 if the growth rate p is sufficiently large 
there exists unique upheaval waves whose speed depends on (21). On the other hand, if p is not 


sufficiently large there is the possibility of retreating waves, where the non-excited state (0,u*(0)) 
invades regions with a high level of rioting activity ( u* , v* (u *)). 

Figure [3] illustrates various numerical experiments of system ([?]). Figure s [3a| and 3b illustrate 
two simulations with the exact parameters (p and (I in region Ha of Figure [lFbut with different 
initial conditions. In particular, the initial conditions used in the simulation illustrated by Figure 
3a is 0(e~ x ), which affect a larger region, leads to a faster upheaval wave than that illustrated 
in Figure 3b In fact, in Figure 3b the initial condition is 0(e~ 3x ). Figure |3c| illustrates another 


simulation with no restriction of information, but with parameters in the bistable region. In this 
case all initial conditions lead to waves of the same speed. Figures |3d]|3f| illustrate simulations with 
some restriction of information. Figure [3d] illustrates a bistable upheaval wave resulting from low 
restriction of information where the excited state invades the non-excited state. As the restriction 
of information is increased the traveling wave slows down and eventually we obtains a stationary 
wave, as is illustrated in Figure [3e[ and even a retreating wave when the restriction of information 
is sufficiently high (Figure [3e[). In the latter case one observes that the non-excited state invades 
the excited state. 
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Level of rioting activity 




(a) a = 0: «o(0) = 5e x (b) a = 0: Mo(0) = 5e 3x 




(c) a = 0, fast transition (d) 0 < a < 1, bistable wave 




(e) 0 < q < 1, stationary wave (f) 0 < a < 1, retreating wave 

Figure 3: Traveling wave solutions for the level of rioting activity with various parameters and initial 
conditions. 


Proof. (Theorem [2]) We rewrite the system given by ( |10a )-( 10b) in vector-matrix form: 

u = F[u\, 


(23) 


where u = (u,v) and F[u\ = First, note that 3>„(?b v) = r'(v)G a (u) > 0 for all 

a G [0,1] since r'(y) > 0. Moreover, 'S> u (u, v) = ( 1+ ^p-i v > 0 for all u,v > 0 and so the system 
( |23| is monotone. Now, the stability of any constant steady-state solution (u c ,v c ) to (23) is given 
by: the sign of the eigenvalues of the matrix: 


F'[(u c ,v c )\ = 


r(v c )G' a (u c ) - 1 r'(v c )G a (u c ) 
—h'(u c )v —( h{u) — kf) 


( 24 ) 


For notational simplicity we denote the trace of (24) by: 


Tr[(u c ,u c )] = r(v c )G' a (u c ) - 1 - (h{u) - k 2 ), 
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and the determinant by 

D[(u c ,u c )] = -( r[v c )G' a {u c ) - 1 ){h(u) - k 2 ) + r'(v c )G a (u c )h(u). 

Case 1: Let us first consider the case when the parameters are chosen from regions Ha, Ilia, or 
Illb of Figure [I] or Ila-IIb of Figure [2j For the non-excited steady-state (0,u*(0)) we have that 

1*1(0, »*(0))] = 1 + e _X.(,)-a) G (.(°) - 1 - (1 - fa) < 0. 

Indeed, when a = 0 we know that G' a ( 0) = 0 and the inequality is clear. For the case when a = 0 
recall that p and /3 are chosen such that 1+e _ |3( ^ (0) _ a) G'q(O) < 2 — k 2 . Similarly, we observe that 

D[( 0 ,u*( 0 ))] = —(r(u*( 0 ))G / Q ,( 0 ) - 1)(1 - k 2 ) + rV(0))G Q (0) > 0 - 

Thus, the eigenvalues of F'[(0, u*(0))] lie on the left-half plane. Now, for the excited state (u*,v* (u*)) 
we have that 

Tr[(u*,v*(u*))] = r(v*(u*))G' a (u*) - 1 - (1 - k 2 ) < 0, 
since G' a (u*) < 0 and by the same token 

= -(r(v*(u*))G' a (u*) - 1)(1 - k 2 ) + r'(v*(u*))G a (u*) > 0. 

Again, from this it is clear that the eigenvalues of F'[(u*,v*(u*))] lie on the left-half plane. For 
parameters in the region Ha on Figure [l] there are only two steady-states. In the remaining regions 
there is a warm state (u w ,v*(u w )) such that 


Tr[(u w ,v*(u w ))] = 


;G' a (u w ) - 1 - (1 - k 2 ) > 0, 


1 + e -P(.v*(u w )-a) 

since G' a (u w ) is sufficiently large so that 1+e -/ 3 ( v P *(u w )-a) G' a (u w ) > 2 — k 2 . Thus, at least one of the 

eigenvalues of F'[(u w ,v*(u w ))] lie on the right-half plane. 

Case 2: Let us now consider the case when the parameters are chosen from region lib in Figure 
[2| In this case, there are two steady-state solutions (0, u*(0)) and (u*,v*(u*)). The analysis for 
the excited state is the same as above so that all of the eigenvalues of the Jacobian evaluated at 
this steady-state lie on the left-half plane. However, for (0,u*(0)) the parameters are such that 
/ 3 (C*(o)— a) ^a( 0 ) > 2 — k 2 so that 

Tr[(0, u*(0))] > 0. 

This implies that at least one of the eigenvalues of the Jacobian evaluated at this steady state lies 
on the right-half plane. 

Thus, the existence and monotonicity of the traveling wave solutions (4>(z), tp(z)) that satisfy 
(20) now follow directly from Theorem 2.1 in f^Sj for 1 and Theorem 2.2 for 2 and we are left to 
verify (i). For this purpose, let U{z ) be a traveling wave solution with speed d satisfying 


U" + c'U' + g a (U) = 0, 

0 < U < u*, 

U{— oo) = u*, U(+ oo) = 0. 


(25) 


Note that such a solution exists and that the speed d depends on F(u*), both in sign and in 
magnitude. We also know that U(z) is monotone decreasing. Now, let V = v*(U) so that 


lim V(z) = v*(u*), lim V(z) = u(0), and 'I ’(V,U ) = 0. 


(26) 
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Without loss of generality assume that c! > 0 and assume for contradiction that c* < 0. Now, 
define w(z, t ) = u(z + ct, t ), b(z, t ) = v(z + ct, t ), and the operator 


( w t - w zz - c*w z - $( 6 , w ), 
{ h - b zz - c*b z - $>(b,w). 


(27) 


For any v := 
(U. V) gives: 


(w,b) which satisfies ( 20 ) it holds that A/V 


0. Applying the operator M to vi = 


f C *U' - g a (U), 

\ -V"-c*V 


(28) 


First, note that since U satisfies (25) we obtain that — U" — g a {U) < 0 and —c*U' < 0. Moreover, 
using (26) we have that — V" — c*V' < 0, which can be seen by integrating the left-hand-side and 
noting that 

—c*(V( oo) — V(—oo)) < 0. 


Thus, we conclude that A/V] < 0, so that (U, V) is a super-solution. To exclude any issues at Too 
let us choose zq € M such that 


U(z + ^o) < 4>{. z ) an d V(z + zq) < i/j(z) for all z < R, 


(29) 


for some R G M and such that U(z + zo) < 6 for some S > 0 and all z > R. In fact, (29) must hold 
for all z E M. To observe this note that 


+g a {4>) < 0, 

U"{z)+g a {U)> 0. 

Thus, if we dehne m(z) = U(z) — cj>(z) then it satisfies m"(z) + q(x)m(z) > 0 where 

, x 9a{U) - g a {4>) 

9(l) = u-4, ■ 

Now, for parameters p and (5 in regions III of Figure [I] or regions II on Figure [2j since U < 5 
we know that q(x) < 0 in the set {m > 0}. An application of the strong maximum principle to 
m"(z) + q(x)m(z ) > 0 on the set {m > 0 } implies that m < 0 as m = 0 on the boundary, which is 
a contradiction, and so we conclude that m < 0 for all z E M. A similar argument works for (f>(z) 
and V. Thus, we know that 

U(z) < <f>(z) and V{z) < ip(z), 

for all z. The assumption that c* < 0 leads to a contradiction, which implies that c* > 0. Similarly, 
we can conclude that if c' = 0 then c* = 0 and if d < 0 then c* < 0 . 

□ 


6 Heterogeneous environments 

As for ecological models, it is clear that the heterogeneities in the environment play a significant 
role on the ignition, spread, and propagation of rioting activity. Indeed, laws, unemployment rates, 
forms of government, etc. vary drastically between countries and even from state to state within 
a country. To take this even further, income inequality is observed between neighborhoods that 
are adjacent to each other. It is therefore important to analyze how these heterogeneities impact 
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the overall failure or success of an upheaval wave. In an effort to gain some perspective in this 
direction we perform various numerical experiments in heterogeneous environments. In general, 
these heterogeneities can be quite complex and difficult to analyze and a good first approximation 
is the so-called patch model, a terminology borrowed from the ecology literature, see for example 
PH EE G3 and references within. A patchy environment refers to a fragmented environment where 
each fragment is relatively homogeneous. Specifically, we choose to work in an environment where 
the heterogeneities are due to variations in the restriction of information of each region. Thus, a(x) 
changes by country or region: we can think of a(x ) as being small for less suppressive states and 
a(x) as being large for states with high restriction of information. Through this simple environment 
we are able to see a possible cause for why riots (or even revolutionary waves) skip certain countries. 
Of course, we cannot conclude that this is the absolute cause; however, it is interesting to note that 
a simple set of assumptions will lead to rioting activities failing in certain regions while it propagates 
globally. 

6.1 Periodically fragmented environments 

Here, we are interested in understanding when a patchy environment leads to the persistence of a 
protest or riot and when it leads to these protests vanishing. For this purpose we study the following 
problem: 


ut = Au + <h(x, u, v), (30a) 

Vt = DAv + ^{u, v), (30b) 


where 3>(x,m, v) is periodic in the spatial variable. For the sake of generality let us assume that 
x G for d > 1 and let L\, L 2 , ■■■, L d > 0 be real numbers. A periodic function g : — > M is a 
function such that g{x 1 , ...,Xk + £&,..., Xd) = g(x 1 ,..., Xd) for any k G {1,..., d} and the period L is 
defined by: 

L = (0, Li) x (0, L 2 ) x • • • x (0, L d ). 

We will assume that $(x,u,v) and ^>(x,u,v) are periodic in x, a natural example being the case 
when <F(x,m, m) = r(v)G a ^(u) — 1 where: 


J Ql X G Si, 

\ a 2 x G L\Si, 


(31) 


where ot\ > 0.2 and S\ C L. Following the notion of stability of the non-excited state in [3], we look 
for the principal eigenvalue, A, of the operator £ = (£ 1 , £ 2 ) defined by: 


£1 := —Am — $> u (x, 0, 0)m — $> v (x, 0, 0)m, 
£2 := -DAv - T u (0, 0)m - \k„(0, 0)m, 


with periodic conditions. Thus, we look for A such that there exists positive functions </>, ip that 
satisfy: 


-A (p - <h u (.x, 0, 0)(p - $ v (x, 0, 0)?/’ = X(p, 

-Aip- ty u (0,0)(p- T,„(0,0)V’ = Xip, (32) 

<t>,fp> 0 periodic, = 1 . 


It can be shown that there exists a unique principal eigenvalue, which is simple, that satisfies (32). 
Indeed, this follows by noticing that replacing A by A + k leads to a cooperative system. Then 
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the strong maximum principle applied and an application of the Krein-Rutman theorem gives the 
result. In the case when A < 0 we say that the non-excited state (0,v*(0)) is unstable and thus we 
expect the existence of a positive steady-state solution. This is stated rigorously in the following 
proposition. 


Proposition 2. If A < 0 there exists a positive and periodic steady-state solution to (30). 


Proof. Since A < 0 we can choose k sufficiently small so that 


and similarly, 


Xk 

&{kip, k(j>) - —cj> > kip$ v (0 ,0) + k(j)$ u (0 , 0), 


Xk 

^{kip, k(j>) - —ip > kip^ v (0 ,0) + kif$ u (p, 0), 


Then multiplying the top equation of (32) by A: gives: 


(p < 0, 

(33) 

■cp < 0. 

(34) 


The inequalities (33) tells us that for k sufficiently small (kip(x), k<p(x)) is a sub-solution. Moreover, 
(u* ,v*(u*)) is a super-solution. Finally, as (u*,v*(u*)) > (kip (rc), kip) then, since the system is 
cooperative, a classical iteration scheme yields a periodic solution (p(x),q(x)) with the property 
that k(p{x) < p{x) < u* and kipfx) < q(x) < v(u*). □ 


As noted earlier the notion of traveling wave solutions is important because these solutions 
represent the spread of high levels of rioting activity and tension into regions with low levels of 
activity. In heterogeneous environments this notion, as introduced in the section [5| does not apply. 
However, the need for such type of solutions still remains true which has lead to the generalized 
notion of such solutions in periodic domains. These solutions are refereed to as pulsating fronts and 
they propagate in a given direction but with a profile that changes periodically [2]. In other words, 
these fronts connect the non-excited state and a periodic steady-state. We observe such solutions in 
numerical simulations and illustrate them in Figure [4j In particular, when there is no restriction of 
information we observe pulsating fronts whose speed depends on the transition parameter. Figure [4a] 
illustrates a fast moving wave that corresponds to a slow transition rate. On the other hand, Figure 
|4b| illustrates a slower moving pulsating wave that corresponds to a fast transition rate. Finally, 
Figure [4c] illustrates a lack of such type of solutions when there is some restriction of information. In 
fact, Proposition [2] provides a hint as to when we expect these pulsating waves to exists. Specifically, 
when the principle eigenvalue satisfying (32) is negative there is hope of a pulsating front existing. 
This is an important open problem which we hope to address in the future. 
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(a) q = 0 and /3 small: m o ( 0) = 5e x 



(b) a = 0 and /3 large 



Figure 4: Pulsating waves in a periodic environment for the level of rioting activity. 

6.2 Gaps in upheaval waves 

Another typical heterogeneous environment studied in various fields, such as ecology, neural net¬ 
works, mathematical criminology, and excitable media (see mmm and references within) is the 
so-called gap (or barrier) problem. In this context, we consider an environment where there are 
three regions each with a different level of restriction of information. For the sake of illustration let 
us imagine that we have three neighboring countries, each state is denoted by Si, for i = 1, 2 ,3 in 
/ C 1. In all of the simulations that we perform the domain of interest is the interval I = [0,15] 
and the middle country, S 2 , has the highest level of restriction of information, that is a = 1. To 
make this precise, we define the heterogeneous restriction of information as follows: 

{ a.\ x G Si, 

1.0 x G 5*2, (35) 

«2 x G S 3 , 

for 0 < ai, 0*2 < 1. Figures [5a|and[5b|illustrate bistable and monostable upheaval waves, respectively, 
which successfully cross the region with high restriction of information. 

Remark 2. This offers a possibility as to why protests or riots fail in certain countries, while being 
globally successful. If the country has a high restriction of information it will not experience protests 
even if the neighboring countries have full force protests. 
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In Figures 5c and 5d we explore the effect of the size of the middle country on whether the 
revolution is able to cross the country with high restriction of information and take off in the 
neighboring country. In turns out that if the country is sufficiently small the upheaval wave will 
successfully cross the country with high restriction of information (still skipping that country), this 
is illustrated in Figure |5d| On the other hand, if the country is too large then this is enough to 
prevent the upheaval wave from going beyond the originating country, this is illustrated in Figure 



(a) Bistable invasive wave 
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(b) Monostable invasive wave 



(c) Bistable blocked wave: large gap (d) Bistable invasive wave: small gap 

Figure 5: Level of rioting activity in a heterogeneous environment with a gap. 


7 Discussion 

This work is a sequel to [5] where we introduce mathematical models for the dynamics of rioting 
activity, which included exogenous and endogenous effects as well as non-geographical connections. 
Here we have presented a way to model restriction of information in a given environment, which 
plays a significant role on whether rioting activity or protests will be successful, in the sense that 
they will spread and last for prolonged periods of time. This model allowed us to explore the effects 
of different heterogeneities in the environment in the persistence/extinction of protests and in the 
existence or non/existence of upheaval waves. 

We conclude that upheaval waves can, in fact, exist even when there is some restriction of 
information. However, when there is no restriction of information these waves can travel arbitrarily 
fast, which is not the case when there is some restriction of information. As the restriction of 
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information increases we also observe the existence of retrieving traveling waves, which essentially 
are solutions with a negative speed. From the point of view of the application we have only begun 
to scratch the surface and there are many questions that remain to be explored. For example, the 
question of whether there are external factors that can help stop certain riots is of interest. 

In addition to vast amount of work that remains to be done from the point of view of the 
application, there is also much mathematical theory that remains to be developed for the systems 
that we introduce here and in |5|. For example, an extremely interesting open problem is to prove 
the existence of traveling wave solutions in the case when the social tension diffuses non-locally. 
Such results hold for single parabolic equitations, under some conditions of the kernel; however, to 
our knowledge the system has not be treated. The existence of pulsating waves in periodic domains, 
which we observe numerically in our simulations, also remains to be proved. A Liouville type result 
for the system of equations is also an interesting avenue of research. 
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8 Appendix 


The appendix is devoted to working through some of the analysis and heuristics used to develop 
Figures [I] and [2] As noted in section 3.2 any constant steady-state solution of ([T]) with k = 0 must 
satisfy: 


v*(u) = 


1 


h(u) — fe 2 

Note that v*(u ) is an increasing function of u such that 

lirn v*(u) = 00 , 


and <J>(u*(tt), u) = 0 . 


i/p 

for u = ( ) — 1. For simplicity, let us define 


H[u ) := &(v*(u),u). 


I. a = 0 case: 


For notational simplicity let us set G{u) = Gq{u). First, we observe that H(u) < pu{ 1 — u) — u 
which is strictly negative for all u > 0 if p < 1. Now, consider the derivative of H(u ) at u = 0 


H'{ 0) 


P 

1 _|_ e -P(v*(0)-a) 


- 1 , 


(36) 
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In particular, for 


P> 


Hp~ i) 

a — u*(0) ’ 


we obtain that H'( 0 ) > 0 . This gives the first curve in Figure |TJ which is defined by /3i(p) = ^ 

and is the curve where H'( 0) = 0. Now, in the region where H'(0) < 0 we need to determine if 
H(u) > 0 for some 0 < u < min {1,12}. For this purpose it is useful to look at the details of H(u) 


H(u ) 


P 


1 _|_ e -/3 (v*(u)-a) 


-1 u- 


P 


1 _|_ e -P(v*(u)-a) 


-ti 



Figure 6: Illustration of the function r{z) with critical tension a = 10 and transition parameter j3 = 1 and 
/3 = 10. 


Note that for v*(u) > a for /3 sufficiently large we have that 1 +e _^(^«( u )_a) ~ p- Thus, for each 
p > 2 there exists a /3 P such that H(u) > 0 for some u < min 1, u with j3 > f3 p . In particular f3 p is 
decreasing. Moreover, if p is sufficiently large then H(u ) > 0 for some u < min ( 1 , u} and all f3. This 
gives the curve 02 (p) in Figure [l] Similarly, we can see that when in the regime where H'( 0) < 0 
but H(u) >0 for some u < min {1,12} so that there are three constant steady-state solutions for p 
in intermediate regions then the quantity F(u*) < 0 as defined in ( 21 ) where u* is the largest zero. 
On the other hand, as p increases F(u*) will become positive. This gives rise to the curve given by 
/?4 (p). Finally, to see where /3s(p) comes from, first notice that from (36) that H'( 0) increases with 
p from zero. On top of that, iT(0) will decrease with /3. 


II. 0 < a < 1 case: 

As soon as the restriction of information parameter is positive we see immediately that -£^( 0 ) < 0 . 
In this case there will generally only be one or three constant steady-state solutions. The possibility 
of two steady-state solutions exists on the curve /3i(p) in Figure [2j The parameters lying on the 
curve f3i(p) leads to a non-zero steady-state solution which is degenerate. In fact, H{u) < 0 for all 
u > 0. For the curve fh(p) we can argue as was done for the a = 0 case. 
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